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In this talk I report my recent study on the shear viscosity of neutron-rich nuclear matter from a 
relaxation time approach. An isospin- and momentum-dependent interaction is used in the study. 
Effects of density, temperature, and isospin asymmetry of nuclear matter on its shear viscosity have 
been discussed. Similar to the symmetry energy, the symmetry shear viscosity is defined and its 
density and temperature dependence are studied. 
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I. INTRODUCTION 



One of the major problems in nuclear physics is to un- 
derstand the properties of nuclear matter under extreme 
conditions. This is related to the basic knowledge of 
the in-medium nuclcon-nuclcon (NN) interaction which 
in the present stage can still hardly be obtained from the 
ab initio theory of the strong interaction, i.e., Quantum 
chromodynamics. Our knowledge on the in-medium NN 
interaction today is mainly developed along two lines. In 
the first line, one starts from the bare NN interaction, 
which has been fitted very well from NN scattering data, 
together with phenomenological three-body interactions, 
so that the in-medium NN interaction and the properties 
of nuclear matter can be obtained through many-body 
theories. In the second line, the starting point is an effec- 
tive in-mcdium NN interaction or Lagrangian, with the 
parameters fitted to the empirical nuclear matter proper- 
ties obtained usually through mean-field approximations. 

Ten years ago, an isospin- and momentum-dependent 
mean-field potential (hereafter 'MDI') was constructed 
to study the dynamics (especially the isospin ef- 
fects) in intermediate-energy heavy-ion collisions to- 
gether with an isospin-dependent Boltzmann-Uehling- 
Uhlenbeck (IBUU) transport model In addition to 
the good description of the empirical nuclear equation 
of states, the momentum dependence of this mean-field 
potential reproduces pretty good the optical potential 
extracted by Hama et al. from elastic proton scatter- 
ing data [2|. The studies using this interaction have 
constrained the nuclear symmetry energy at both sub- 
saturation and suprasaturation densities In ad- 
dition to the dynamics of heavy-ion collisions, the MDI 
model has also been used to the study the thermodynam- 
ical properties of nuclear matter @, 0] ■ It was recently 
found that the isospin- and momentum-dependent poten- 
tial can be derived from an effective interaction with a 
density-dependent two-body interaction and a Yukawa- 
type finite-range interaction using Hartree-Fock calcula- 



tion [8[ . The MDI model thus serves as a useful effective 
in-medium interaction. 

In the past few years, the shear viscosity of the quark- 
gluon plasma (QGP) formed in relativistic heavy-ion col- 
lisions has attracted special attentions. From the study 
with a viscous hydrodynamical model it was found 
that the strong-interacting QGP behaves like a nearly 
ideal fluid, i.e., its specific shear viscosity is only a little 
larger than the KSS boundary [l^]. Up to now large ef- 
forts have been devoted to study the shear viscosity of 
QGP pl] - [l^ and hadron resonance gas [isl - fisj formed in 
relativistic heavy-ion collisions, while there are only a few 
studies on the shear viscosity of nuclear matter formed 
in intermediate-energy heavy- ion collisions [T9l - [23| . Even 
few studies are related to the isospin effects on the shear 
viscosity of nuclear matter 12411. In the present talk I 
will discuss my recent study [25| on the shear viscosity 
of nuclear matter using the MDI model mentioned above 
from a relaxation time approach, which gives an intuitive 
picture how the shear viscosity changes with the density, 
temperature, and isospin asymmetry of nuclear matter. 



II. SHEAR VISCOSITY FROM A RELAXATION 
TIME APPROACH 



The system concerned here is an isospin asymmetric 
nuclear matter with uniform neutron and proton density 
Pn and pp, respectively, and the nucleons are thermalized 
with temperature T. The flow flcld u is static in the z 
direction and its magnitude is linear in the coordinate 
X, i.e., Uz = cx and = Uy = 0. In the rest frame 
nucleons move with the flow flcld and follow Fermi-Dirac 
distribution n* in the equilibrium state. In the lab frame 
the equilibrium distribution is a simple boost by the flow 
field compared with that in the rest frame, denoted as nP . 
Due to NN collisions, the real distribution may be slightly 
away from the equilibrium distribution and is denoted as 
n, and the deviation from the equilibrium distribution 



5n ~ 



is much smaller than 
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The shear force between flow layers per unit area by 
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definition can be written as 



F 
A 



< {p^ - mUz)prVx > 



(1) 



the momentum distribution n,- = n° + Srir to calculate 
the average and taking into account that the equilibrium 
momentum distribution n° is even in px, Eq. ([1]) can be 
fmthcr written as 



In the above, r = n or p denotes the isospin degree of 
freedom, prVx is the number of nucleons moving between 
layers per unit time per unit area, and Pz—muz is the mo- 
mentum transfer per nucleon in the z direction. The nu- 
cleon velocity in the x direction Vx can be further written 
as I'x ~ Px/m*, with m* being the effective mass. Using 



J 



-^l^dJ{pz-muz)—Sn.^, 



(2) 



where d = 2 is the spin degeneracy. 



In the following I will calculate Srir by linearizing the isospin-dependent BUU equation as follows 



dt 



V-VrMPl)-yrUr-V,nr{pi) = -{d--) J ^^^jT,,, 

d^P2 d^p'i d^p'2 



X [nr{pi)n^r{P2)(l - nr{p[)){l - n-r{p'2)) - ?v(p'l)'l-r(P2)(l - '^r(Pl))(l 

X {2nf5^'Hpi+p2-I^^-I^2)- 



-iP2))] 



(3) 



In the above, T is the transition matrix, the degeneracy d— 1/2 takes the double counting of identical nucleon collisions 
into consideration, and 1 — n is from the Pauli blocking effect. Replacing n with nP in the first-order approximation, 
the left-hand side can be expressed as 



dnrjpi) 
dt 



duz PzPx dn° 
dx p dp 



(4) 



by using the properties of n*'. Keeping only the Snr{pi) term, the right-hand side of Eq. 1^ can be expressed as 
6nr{pi)/Tr{pi), where Tr{pi) is the relaxation time, i.e., the average time between two collisions for a nucleon with 
isospin T and momentum pi , and it can be written as 



Tr(Pl) 



^(Pl) 



diff 



iPl)' 



(5) 



where Tt 



sarne{diff) 



(pi) is the average time for a nucleon with momentum pi to collide with other nucleons of same 



(different) isospin, and they can be calculated respectively from 
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p\dp2d cos 6*12^ cos 9 
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dn 
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Tt 



dtff 
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(7) 



In the above 9i2 is the angel between pi and p2, and 9 
is the scattering angel between the total momentum and 
the relative momentum of the final state. In free space 
the pp and np scattering cross sections are isotropic and 



they can be respectively parameterized as (26| 

(Jppinn) = 13.73 - 15.04/w + 8.76/v^ + 68.67w^, (8) 
cr„p = -70.67 - 18.18/w + 25. 26/^2 + 113. 85w, (9) 

where the cross sections are in mb and v is the velocity 
of the projectile nucleon with respect to the fixed target 



3 



nucleon. It is worth to note that in the most probable 
colhsion energies the np scattering cross section is about 
three times the pp scattering cross section. In nuclear 
matter, the in-medium NN scattering cross sections are 
modified by the in-medium effective mass in the form of 



medium 
'NN 



(^NN 



^^NN 
fJ-NN 



(10) 



where /inn (/^jvat) is the free-space (in- medium) reduced 
mass of colliding nucleons. 

Once the relaxation time Tr{p) is known, SriT^p) can 
be calculated from 



5nr{p) = Tr{p) 



duz PzPx dn° 
dx p dp 



(11) 



Using the definition F/A = ~ij{duz/dx), the shear vis- 
cosity can be calculated from Eqs. ([2]) and (fTTj) in terms 
of the local momentum distribution nt as 



Trip) 



plpldn* (pp 
pm* dp (2-kY ' 



(12) 



by setting the magnitude of the velocity field to be in- 
finitely small. Note that from Eq. ([T2|) the shear viscosity 
is related to the local momentum distribution near the 
Fermi surface. 



III. RESULTS AND DISCUSSIONS 

Figure [1] displays the density, temperature, and isospin 
dependence of the relaxation time. In neutron-rich nu- 
clear matter, t!^^^^ is larger while r^"™'^ is smaller com- 
pared to that in symmetric nuclear matter as a result 
of less frequent np collisions and more frequent nn colli- 
sions. For the similar reason, r*'""'^ is larger while Tp*-'"'^ 
is smaller compared to that in symmetric nuclear mat- 
ter. From Eq. ([5]), the total relaxation time is determined 
by r'^'^ff which is always smaller than r^"™*^ due to the 
larger np cross section than pp (nn) cross section in the 
most probably collision energies. Thus, neutrons have a 
larger relaxation time than protons in neutron-rich nu- 
clear matter. It is seen in Panel (d) that the relaxation 
time decreases with increasing temperature due to more 
frequent collisions at higher temperatures. In addition, 
at lower temperatures the relaxation time peaks around 
the Fermi momentum, indicating a strong Pauli blocking 
effect for nucleons near the Fermi surface. 

Results of the shear viscosity rj and specific shear vis- 
cosity rj/ s, where s is the entropy density, are shown in 
Fig. [5] The temperature dependence of the shear vis- 
cosity is similar to that in Ref. [l^] at different densi- 
ties, while rj increases with increasing density especially 
at lower temperatures due to the strong Pauli blocking 
effect. The specific shear viscosity decreases with in- 
creasing temperature, and it is similar in both magni- 
tude and trend to those obtained from BUU calculations 
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FIG. 1: (color online) Panel (a), (b), (c): Relaxation time for 
neutrons and protons as a function of nucleon momentum in 
symmetric (5 = 0) and asymmetric (5 = 0.5) nuclear matter 
at saturation density and temperature T = 50 MeV; Panel 
(d): Relaxation time as a function of nucleon momentum in 
symmetric nuclear matter at different densities and tempera- 
tures. 
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FIG. 2: (color online) Density and temperature dependence 
of the shear viscosity ((a), (b)) and specific shear viscosity 
((c), (d)) for symmetric (5 = 0) and asymmetric (5 = 0.5) 
nuclear matter. 



using the Green-Kubo formula [23l |. It is interesting to 
see that at higher temperatures the specific shear viscos- 
ity is about 4 ~ 5 times the lower limit from Ads/CFT 
calculation [l3|, which is already close to that of QGP ex- 
tracted from the study using a viscous hydrodynamical 
model 0. At lower temperatures the specific viscosity 
increases with increasing density due to the Pauli block- 
ing effect, while at higher temperatures the dependence 
on the density is rather weak. 

Due to the sharper momentum distribution of neu- 
trons compared to that of protons in neutron-rich nu- 
clear matter, the total shear viscosity is dominated by 
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FIG. 3: (color online) Density and temperature dependence of 
symmetry shear viscosity ((a) and (b)) and symmetry specific 
shear viscosity ((c) and (d)). 

neutrons which have a longer relaxation time in asym- 
metric nuclear matter compared to that in symmetric 
nuclear matter. This is confirmed in Fig. [2] that both 
the shear viscosity and the specific shear viscosity are 
larger in neutron-rich nuclear matter. In addition, it 
was seen [2^ that both the shear viscosity and specific 
shear viscosity satisfy the parabolic approximation with 
respect to the isospin asymmetry, i.e., 

v{p,T,5) « 7^{p,T,5 = Q)+ii,y,^{pT)5\ (13) 
CI) (p, T, 5) « (^) (p, T,6 = Q) + (^) (p r)4i4) 

Vs/ Vs/ Vs / sym 



Similar to the symmetry energy, the second-order coeffi- 
cient can thus be defined as the symmetry shear viscos- 
ity or the symmetry specific shear viscosity. The den- 
sity and temperature dependence of them are shown in 
Fig. [3l It is seen that both the symmetry shear viscos- 
ity and symmetry specific shear viscosity decrease with 
increasing temperature. At lower temperatures, both of 
them increase with increasing density. At higher temper- 
atures, the density dependence is rather weak, yy^ym and 
(f )sym ^''^ important quantities in understanding trans- 
port properties of neutron-rich nuclear matter, and they 
deserve further studies in the future. 



IV. SUMMARY AND OUTLOOK 

Using a relaxation time approach, I studied the shear 
viscosity and specific shear viscosity of hot neutron-rich 
nuclear matter as that formed in intermediate-energy 
heavy-ion collisions by using an isospin- and momentum- 
dependent interaction. It is found that the specific shear 
viscosity decreases with increasing temperature, and it 
increases with increasing density at lower temperatures 
due to the strong Pauli blocking effect. Furthermore, 
both the shear viscosity and specific shear viscosity are 
found to increase with increasing isospin asymmetry of 
nuclear matter and roughly satisfy the parabolic approx- 
imation. The second-order coefficient in the expansion of 
the isospin asymmetry, which is defined as the symmetry 
shear viscosity or the symmetry specific shear viscosity, 
has also been studied. 



[1] C.B. Das, S. Das Gupta, C. Gale, and B. A. Li, Phys. 

Rev. C 67, 034611 (2003). 
[2] S. Hama, B.C. Clark, E.D. Cooper, H.S. Sherif, and R.L. 

Mercer, Phys. Rev. C 41, 2737 (1990). 
[3] L.W. Chen, CM. Ko, and B.A. Li, Phys. Rev. Lett. 94, 

032701 (2005). 

[4] B.A. Li and L.W. Chen, Phys. Rev. C 72, 064611 (2005). 
[5] Z.G. Xiao, B.A. Li, L.W. Chen, G.C. Yong, and M. 

Zhang, Phys. Rev. Lett. 102, 062502 (2009). 
[6] J. Xu, L.W. Chen, B A. Li, and H.R. Ma, Phys. Lett. 

B650, 348 (2007). 
[7] J. Xu, L.W. Chen, B. A. Li, and H.R. Ma, Phys. Rev. C 

77, 014302 (2008). 
[8] J. Xu and CM. Ko, Phys. Rev. C 82, 044311 (2010). 
[9] H.C Song, S.A. Bass, U. Heinz, T. Hirano, and C Shen, 

Phys. Rev. Lett. 106, 192301 (2011). 
[10] P.K. Kovtun, D.T. Son, and A.O. Starinets, Phys. Rev. 

Lett. 94, 111601 (2005). 
[11] A. Peshier and W. Gassing, Phys. Rev. Lett. 94, 172301 

(2005). 

[12] A. Majumder, B. Miiller, and X.N. Wang, Phys. Rev. 

Lett. 99, 192301 (2007). 
[13] Z. Xu and C. Greiner, Phys. Rev. Lett. 100, 172301 



(2008) . 

[14] J.W. Chen, H. Dong, K. Ohnishi, and Q. Wang, Phys. 

Lett. B685, 277 (2010). 
[15] A. Muronga, Phys. Rev. C 69, 044901 (2004). 
[16] J.W. Chen and E. Nakano, Phys. Lett. B647, 371 (2007). 
[17] N. Demh and S.A. Bass, Phys. Rev. Lett. 102, 172302 

(2009) . 

[18] S. Pal, Phys. Lett. B684, 211 (2010). 
[19] P. Danielewicz, Phys. Lett. B146, 168 (1984). 
[20] L. Shi and P. Danielewicz, Phys. Rev. C 68, 064604 
(2003). 

[21] J.W. Chen, Y.H. Li, Y.F. Liu, and E. Nakano, Phys. Rev. 

D 76, 114011 (2007). 
[22] S. Pal, Phys. Rev. C 81, 051601(R) (2010). 
[23] S.X. Li, D.Q. Fang, Y.G. Ma, and CL. Zhou, Phys. Rev. 

C 84, 024607 (2011). 
[24] H.F. Zhang, U. Lombardo, and W. Zuo, Phys. Rev. C 

82, 015805 (2010). 
[25] J. Xu, Phys. Rev. C 84, 064603 (2011). 
[26] S.K. Charagi and S.K. Gupta, Phys. Rev. C 41, 1610 

(1990). 



